The development of pair photon cascades initiated by high energy electrons above a pulsar polar cap is simulated numerically. The calculation uses the energy of the primary electron, the magnetic field strength, and the period of rotation as parameters and follows the curvature radiation emitted by the primary, the conversion of this radiation e(+) - e(-) pairs in the intense fields, and the quantized synchrotron radiation by the secondary pairs. A recursive technique allows the tracing of an indefinite number of generations using a Monte Carlo method. Gamma ray and pair spectra are calculated for cascades in different parts of the polar cap and with different acceleration models. It is found that synchrotron radiation from secondary pairs makes an important contribution to the gamma ray spectrum above 25 MeV, and that the final gamma ray and pair spectra are insensitive to the height of the accelerating region, as long as the acceleration of the primary electrons is not limited by radiation reaction.